Os 



O 



o 



The Shape of Inflated Vesicles 

G. Gompper^ and D.M. Kroll^'^ 

( a ) Sektion Physik der Ludwig-Maximilians-Universitat Miinchen 
Theresienstr. 37, 8000 Miinchen 2, Germany 

W AHPCRC, University of Minnesota 



w , 1100 Washington Avenue South 

00 ' 

(N) ! Minneapolis, MN 55415, USA 

> 

g i Abstract 



The conformation and scaling properties of self-avoiding fluid vesicles with zero 



q^ • extrinsic bending rigidity subject to an internal pressure increment Ap > are 

studied using Monte Carlo methods and scaling arguments. With increasing pres- 
sure, there is a first-order transition from a collapsed branched polymer phase to 
an extended inflated phase. The scaling behavior of the radius of gyration, the 
asphericities, and several other quantities characterizing the average shape of a 

x 

vesicle are studied in detail. In the inflated phase, continuously variable fractal 
shapes are found to be controlled by the scaling variable x = ApN 3l/ ' 2 (or equiva- 
lently, y = < V >/N 3u ' 2 ), where N is the number of monomers in the vesicle and 
V the enclosed volume. The scaling behavior in the inflated phase is described by 
a new exponent v = 0.787 ± 0.02. 
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1. Introduction 

In the absence of lateral tension, the shape of a fluid membrane is determined 
by its elastic bending energy [1,2]. At length scales small compared to the per- 
sistence length [3] £ p = aoexp(4rft//csT), where ao is a microscopic length of the 
order of the monomer size and k is the bending rigidity, equilibrium configurations 
minimize the elastic bending energy, and thermal fluctuations can be ignored. At 
longer length scales, however, fluctuations lead to a decrease in the effective bend- 
ing rigidity [4] so that membranes of linear dimension L ^> £ p are expected to 
behave as crumpled objects characterized by the absence of long-range orienta- 
tional order. At these length scales, where the effective rigidity is negligibly small, 
there is no energy scale and fluctuations determine the equilibrium shape. 

Closed membranes or vesicles, usually of spherical topology, are of particular 
interest. For vesicles of radius R <C £ p , it has been shown [5-7] that the intricate 
shapes observed in experiment correspond to equilibrium configurations which 
minimize the Helfrich elastic bending energy [5] 

H = 2k I dS (H - H ) 2 (1) 

subject to the constraints of constant area and enclosed volume, where H = 
^(1/ R\ + I/-R2) is the mean curvature (expressed in terms of the principle radii of 
curvature R\, R%), and H is the spontaneous curvature. Here, we are interested 
in the limit of small bending rigidities [8], where £ p <C R. This case has been 
studied previously in considerable detail in (d = 2)-spatial dimensions, where the 
"vesicle" is a closed ring polymer [9-12]. The shape of these polymer rings can be 
characterized by the thermal averages of the eigenvalues A2 > Ai of the moment 
of inertia tensor. In particular, it has been shown [10] that the value of the shape 
parameter E =< Ai > / < A2 > changes continuously with the pressure increment 

Ap = p in ~ p out . (2) 

For Ap = 0, the polymer ring is a self-avoiding random walk, with E ~ 0.39. For 



Ap 7^ 0, E is a universal function of the scaling variable 

x = pN 2u , with p = Apa 2 /kBT, (3) 

where N is the number of monomers in the ring, a is a length scale of the order of 
the interparticle spacing, and v = 3/4 is the self-avoiding random walk exponent 
in two dimensions which describes the dependence of the radius of gyration, R g , 
on the number of monomers via 

< R 2 g > ~ N 2u . (4) 

For Ap > 0, the rings are inflated and approach a circular shape for large Ap, so 
that E(x) — ► E + = 1 for x — ► oo. For Ap < 0, the rings collapse to a branched 
polymer shape, with E(x) — ► E~ ~ 0.23 in the limit a; — ► — oo. 

In this paper we consider the corresponding problem for self-avoiding fluid 
vesicles in d = 3 spatial dimensions. Our analysis is based on an extensive Monte 
Carlo study of a simple string-and-bead model for randomly triangulated vesicles. 
The surface is modeled by a triangular network of N hard-sphere particles of diam- 
eter a = 1 connected by tethers of maximum extension £q = \[2 [13]. No explicit 
curvature term is included in the model. The fluid character of the membrane is 
incorporated by dynamically updating the triangulations. Details concerning the 
simulation procedure can be found in Refs. 14-16. 

Using this model we have recently shown [16-19] that the behavior of fluid vesi- 
cles in d = 3 spatial dimensions is drastically different from that of ring polymers 
in two-dimensions. For example, already at Ap = a vesicle exhibits a branched 
polymer conformation [16,17,20-23] characterized by the scaling laws [24] 



^> ~ JV "'", 

(5) 



< Rg > ~ N^ 



< V > ~iV^", 

with Ub p = 1. As Ap increases, there is a slow inflation until at some critical 
pressure difference Ap*(N) > 0, there is a first-order transition to a high pressure 



inflated phase [18]. Increasing the pressure still further causes the vesicle to slowly 
inflate until it becomes, ultimately, spherical. In this paper, we characterize the 
inflated phase in some detail. In Sec. 2, scaling laws for the radius of gyration, 
the enclosed volume, and the width of the vesicle wall are derived using a gener- 
alization of de Gennes' blob picture. The scaling behavior of the critical pressure 
difference Ap*(N) at the first-order transition is discussed in Sec. 3, and data for 
the radius of gyration and enclosed volume in the inflated phase are analyzed. A 
simple argument is given which indicates that Ap* ~ iV -1 ' 2 . Sec. 4 contains a 
discussion of the behavior in the two-phase region, and data taken in both the 
constant pressure and constant volume ensembles are compared. A scaling anal- 
ysis of various quantities characterizing the vesicle shape in the inflated phase is 
presented in Sec. 5; the paper closes with a brief summary and discussion. 

2. Scaling 

The scaling behavior in the inflated phase can be derived using a simple 
generalization of the blob picture of de Gennes [11,25-27,18]. Consider a (d — 1)- 
dimensional crumpled membrane embedded in <i-dimensional space. In the blob 
picture, one envisages the membrane breaking up into a set of Nj, weakly stretched 
blobs of area proportional to £^ _1 = hsT/a, where £ CT is the tensile length, under 
the influence of a uniform extensional tension a applied to the membrane perime- 
ter. Since the blobs themselves are only weakly stretched, it is further assumed 
that 

St X ~ M£, (6) 

where Mb = N/Nf, is the number of monomers in a blob. It is then argued that 
the blobs become independent on length scales much larger than £ CT so that the 
total projected surface area < Af, > is proportional to ^~ 1 (N/Mf ) ) ~ No 1 /" -1 . 
If an inflated vesicle is regarded as a (hyper-) spherical bubble of radius R with 
a surface tension 7, one has Ap — 2 r y/R. Identifying the surface tension with the 



stretching tension a, and taking < A\, >~ R d 1 , we arrive at the scaling form 

R g ~ 2W+ (7) 

for the radius of gyration of an inflated vesicle, with p = Apa d /ksT, and the 

exponents 

l-i/ 



A, - 1 (g) 



^+ 



v is a new exponent which characterizes the scaling behavior within a blob in the 
inflated phase. 

Eq. (7) is the large extension limit of a more general scaling form. Assuming 
that the radius of gyration R g and the tensile strength £ CT are the only two relevant 
length scales in the problem, it is natural to expect the projected area At to scale 
as 

A b = R^Shfa/Rg) = N v SI(<tN v ), (9) 

where the last relation follows from the definition R d ~ x ~ N u , and £^ -1 ~ a -1 . 
For large a, A b should vary linearly with N, so that 

Q(z) ~ z~~ l for z — >• oo. (10) 

This implies A^ ~ Na^ 1 ~ u '' 1 ' in this limit. For an inflated vesicle we have (see 
discussion preceding Eq.(7)) a ~ pR g so that 

A b ~R*- 1 ~N(pR g )( 1 -''V v (11) 

in the large extension limit. Solving this equation for R g , we arrive at (7) with the 
exponents (8). From (9), we see that the full scaling form in the inflated phase is 

< R 2 g >= p 2u N 2v +S p (pN dv ^ d - 1 ^). (12) 
Similarly, the mean volume can be shown to scale as 

< V >= p duj N du +T p {pN dv ^ d -^). (13) 
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Finally, Eq. (13) suggests that the scaling behavior of the mean square radius 
of gyration of vesicles fluctuating at fixed enclosed volume (instead of at fixed 
pressure) scales as 

< R 2 g >= V 2/d E v (VN- du/(d - 1) ). (14) 

in the single phase inflated regime. Note that for d = 2, these results are identical 
with the those of Ref. 11 for polymer rings. 

We can go further and calculate the radial extension of the highly convoluted 
surface of the vesicle. The radial thickness £ should scale like the projected area, 

Eq.(9), 

f*" 1 = N v Q{aN v ). (15) 

In the large tension limit, £ should become independent of N, so that 

Q(z) ~ z~ l for z — > oo. (16) 

This implies that 

£ ~ (pRgyw-v ~ p -- +N --+/(d-i) ( 1? ) 

in this limit. The thickness £ described by (17) is just the radius (6) of the blobs in 
the blob picture. It should therefore be considered to be the intrinsic thickness of 
the vesicle walls. In particular, it does not take into account the effect of capillary 
waves on length scales larger than £. These excitations can be easily estimated 
using an effective Hamiltonian for a nearly flat interface of projected area A^ (with 
periodic boundary conditions). We introduce a coarse grained variable u(r) which 
measures the distance of the membrane from some reference plane. On length 
scales larger than £, the lateral tensions dominates the fluctuations, so that 

We// = | / d d -\{Vu)\ (18) 

1 JA h 

with an implicit short length scale cutoff £. The amplitude of the fluctuations is 



calculated with the help of the equipartition theorem, as usual, with the result 



< u 2 > 



d d ~ 1 q k B T 

d-l rrnl 



RJ 1 <*<*->■ ( 27f ) d_1 a( i 



^R g ~ d for d < 3 

hfl\n(R g /Z) forrf = 3 



(19) 



~ < ■ 9 



Inserting the results for a and R g , we arrive at 

(pJV^/(d-i))-V(^-i) forrf<3 



< u 2 > 



<Rj> ~ I (piy 3 -/ 2 )" 1/(3l/_1) Jin (pJV 3 "/ 2 ) f or d = 3 



(20) 



so that capillary wave broadening dominates the thickness of the vesicle walls for 
all d < 3. However, for d = 3, which is the case of interest, the effect of capillary 
waves is very weak, and can only be observed for very large vesicles. 

3. Radius of Gyration and Volume 

The behavior of the average volume < V > of a vesicle as a function of Ap is 
shown in Fig. 1 for three different system sizes, iV = 47, 127, and 247. Each data 
point corresponds to an equilibrium average over 20 — 40 million Monte Carlo steps 
per monomer (MCS). The evolution of a discontinuity in < V > with increasing 
N is clearly visible. Since the transition pressure is difficult to determine from 
< V(p) > itself, we have also calculated the probability distribution P(V). For 
iV = 47, the distribution has a single sharp peak for low and high pressures, but 
broadens and exhibits a weakly bimodal structure near Ap* = 1.05 ± 0.05, which 
we identify with the location of the transition. The distribution for iV = 127 
is shown in Fig. 2 for various pressures near the transition. For Ap = 0.50, a 
bimodal distribution is found, with a somewhat larger peak at large volumes. We 
estimate that the transition occurs at Ap* = 0.47 ±0.02 in this case. For iV = 247 
it is more difficult to determine the location of the transition, since the vesicle 
jumps back and forth between the inflated and the branched polymer states so 
rarely that the equilibrium probability distribution for V could not be determined 
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with sufficient accuracy. However, by studying the time dependence of V for 
Ap = 0.30 and Ap = 0.35, see Fig. 3, we conclude that the transition occurs at 
Ap* = 0.34 ± 0.01 for N = 247. From the three vesicles sizes we have studied, it is 
difficult to determine the iV-dependence of Ap* . Nevertheless, we present a simple 
argument at the end of this Section which suggests that Ap* ~ N~ x ' 2 . This is in 
excellent agreement with our estimates for the transition pressures quoted above 
for N = 127 and N = 247. 

As already shown in Ref. 18, the data for < V > in the inflated phase 
(Ap > Ap*) scale according to (13) with an exponent v = 0.787 ± 0.020. In fact, 
all our data for both the volume and the mean square radius of gyration in the 
inflated regime scale for this value of v. The scaling function F for the volume, 

F(pN 3iy / 2 ) =<V> /(p 3u N 3v +), (21) 

and the scaling function £y(< V > /N 3u ' 2 ) for the mean square radius of 
gyration, see Eq. (14), are shown in Fig. 4. At very large inflations, scaling breaks 
down. The deviation of the last (large Ap) data point for each vesicle size in Fig. 
4a is due to this effect and thus delimits the scaling region. For large scaling 
argument, F(y) must approach a constant value. It is interesting to note that this 
occurs rather rapidly, so that significant deviations can be observed only for small 
pressures and vesicle sizes. In fact, for N — ► oo, this asymptotic limit is is attained 
in all of the inflated phase, since as argued below, Ap* scales as N~ x ' 2 for large 
N so that ApN 3v / 2 -> oo for all Ap > Ap* . 

In Fig. 4b, as well as in Fig. 7 below, we have replaced N by (N — No) in 
the scaling function in order to hasten convergence when N — »■ oo. This parame- 
terization incorporates the leading corrections to scaling. 

Given these results, there is a simple argument which indicates that the tran- 
sition between the collapsed and inflated phases occurs at a critical pressure dif- 
ference Ap* which scales as 

Ap* ~ iV" 1 / 2 , (22) 



independent of the value off. For Ap = 0, the vesicle exhibits a branched polymer 
conformation. To leading order, the configurational entropy S(N) of a branched 
polymer of length N is given S(N) = Nln(z) — 9\n(N), where z is a nonuniversal 
quantity independent of N and 9 — 3/2 [24] . The free energy in the deflated phase 
therefore scales as 

^T- = -N ln(z) - cApN (23) 

to leading order in N and Ap, where c is a constant of order unity. In writing (23) 
we have included the leading pressure dependent contribution [18] Ap < V >~ 
ApN. In the inflated phase, on the other hand, we know that < V >= —dpj-'/dAp 
scales as V + Ap 3uj N 3u + [28]. The leading contribution to the free energy in the 
inflated phase is, therefore, 

V + 

f3F+ = -^(Apf^N 3 ^. (24) 

3u + 1 

The transition occurs when the free energies of the deflated and inflated phases 
become equal, i.e. for 

iV~ (Ap*) Suj+1 N 3iy +, (25) 

or, equivalently, Ap* ~ iV -1 / 2 , independent of v. 

4. Two-Phase Region 

We have performed simulations in both the constant-pressure and 
constant-volume ensembles. Our motivation is twofold. First, we want to check 
the equivalence of the two ensembles for large vesicles. This is important because 
experiments are usually performed at constant volume [7,29]. Second, we want to 
see what a vesicle looks like in the 'two-phase' region. Namely, is there two-phase 
coexistence between a branched polymer region and an inflated part of the vesicle? 

A few typical configurations for iV = 247 and V = 120 are shown in Fig. 5. 
They demonstrate that the vesicle does indeed exhibit two-phase coexistence, with 
a tree-like part attached to a quasi-spherical region. Also, we find that iV = 247 is 
large enough to give essentially identical results in the inflated phase for the two 
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ensembles. This can be seen in Fig. 4b, where the scaling function of the mean 
square radius of gyration is plotted for a range of volumes covering the inflated 
phase, the two-phase region, as well as the collapsed, branched polymer phase. 
Data obtained in both ensembles are included in this figure. Note that the data 
appear to scale over the whole range of volumes. We will see this remarkable 
behavior for other scaling functions in Section 5. We want to emphasize, however, 
that no true scaling can occur in the two-phase region. 

For fully extended spherical vesicles, the volume and the radius of gyration 
are related by V = ^Rg- Thus, 

n = - < V >< R 2 g >" 3/2 (26) 

is a measure for the deviation of the shape of the vesicle from an ideal sphere. 
Our results for 7r(< V > N~ 3u / 2 ) are shown in Fig. 6. For V large at fixed N, 
tv — ► n, as expected. Fig. 6 shows again that the data appear to scale not only in 
the inflated phase, but in the two-phase region as well. 

5. Mean Vesicle Shapes in Inflated Phase 

In order to characterize the shape of the vesicles, we have studied the eigen- 
values Ai < A2 < A3 of the moment of inertia tensor, 

^=^j£ff[>?-r7][rf->f], (27) 

where rj is the position of monomer i, and qi is its coordination number. In 
particular, we have analysed the asphericities, Y\ =< A1/A3 >, T2 =< A2/A3 >, 
as well as the quantities [17,30,31] 

d < Tr T 2 > 
Ad= d-K(TrTr> (28) 



and 



d 2 < Tr T 3 > 

Sd= (rf-l)(d-2)<(TrT)3>' (29) 
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where T a p = T a p — A5 a/ g, and A = h X^i=i -V — ^d < 1 is a normalized measure 

of the anisotropy. The sign of — ^1^3 < Sd < 1 determines whether the object is 

oblate (Sd < 0), or prolate (Sd > 0); its magnitude is a measure of the strength 

of the anisotropy. For d = 3, one has 

A = < Aj + Aj + A| - (AiAg + A 2 A 3 + AgAQ > 
3 <(A!+A 2 + A 3 ) 2 > 

and 



(30) 



<(A 1 -A)(A 2 -A)(A 3 -A)> 

53 = 2~^A^ • (31) 

Let us first consider the averaged eigenvalues < Aj >, i = 1,2,3, of the 

moment of intertia tensor. In the inflated phase, the deviation of the vesicle from 

a completely spherical shape is determined by the fluctuations studied in Section 

2. To leading order, we therefore expect 



<A,>^- 



1 r ' 



< R 2 > + -CiV^ " i2 



<u A > 



(32) 



so that [32] 



" <A ' "' (y) = i + ciy-w-w, (33) 



< Ai + A 2 + A 3 > 
with y =< V > N~ Sl/ / 2 . Since £\ < Xi >=< R 2 >, one has ^. a = 0. With the 

critical exponent v = 0.787 determined above, we therefore expect [3 < A^ > / < 

R > —1] to decay asymptotically as y -1 - 565 . The Monte Carlo data are shown 

in Fig. 7. The eigenvalues < Ai > and < A 3 > exhibit the expected asymptotic 

power law behavior with an exponent which is very close to the value quoted 

above. [1 — 3 < A 2 > / < R 2 >] , on the other hand, decays much faster for large 

y, approximately as y~ 305 . Indeed, the data are consistent with c 2 ~ 0. Since we 

expect, in general, that corrections to spherical shape in the large inflation limit can 



be expressed as a power series in J < u 2 > / < R 2 >, the asymptotic behavior in 
this case should be given by the next to leading term, which is y~ 2 /( 3 ~ 3 ") = y -313 , 
in good agreement with our numerical result. From Fig. 7, the constants q in Eq. 
(33) are found to be 

ci = -0.012, c 2 = 0, and c 3 = 0.013, (34) 
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consistent with the requirement £V Cj = 0. 

The asphericities can be analysed in a similar way. Under the assumption 
that T lj2 =< Ai j2 /A 3 >~< Ai ;2 > / < A 3 >= E i)2 , we obtain [32] 

l-r 1 , 2 (2/) = -(c 1 , 2 -c 3 )y- 1 /(3-3-) (35) 

with y =< V > N~ 3u ' 2 . It has been shown in Ref. 17 that for Ap = the I\ and 
Sj differ by very little, so that this assumption seems justified. The data shown in 
Fig. 8 are in reasonable agreement with the scaling form (35), with an amplitude 
which is indeed close to the difference of the Cj's. Note, however, that Eq. (35) is 
only an approximation, and that in general, one would expect 1 — 1\ to scale with 
an independent exponent. In fact, it appears that the asymptotic slope of 1 — I\ 
does differ somewhat from (35). 

Finally, the shape parameters A 3 and S3 are shown in Fig. 9. The data again 
fall onto a single curve when the scaling variable y =< V > N~ 3u ' 2 is used. The 
asymptotic behavior in this case is found to be 

A 3 (y)~3T 7A 

(36) 

#3(2/) ~?/ 7S , 

with 7a — 3.31, and 75 ~ 6.22. Since these quantities are also averages of the 
eigenvalues A^, it is tempting to conjecture that 

7A ~2/(3-3z/) and 7s ^ 4/(3 - 3z/). (37) 

Note that S3 > 0, but very small, so that the vesicles are prolate on average. 

6. Summary and Discussion 

The scaling behavior of fluid vesicles subject to an internal pressure incre- 
ment Ap > has been studied. Clear evidence has been presented for a phase 
transition from a low-pressure branched-polymer phase to a high-pressure inflated 
phase. Several quantities characterizing the vesicle shape, such as the average 
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eigenvalues of the moment of intertia tensor and the asphericities, have been de- 
termined using Monte Carlo methods. In the inflated phase, all quantities are 
found to be universal functions of the scaling variable x = ApN 3u ' 2 (or equiva- 
lently, y = < V >iV~ 3l/ / 2 ), with v = 0.787±0.02. The asymptotic scaling behavior 
for large inflation can be understood using simple scaling arguments. 

Much still needs to be learned regarding the properties of the inflated phase, 
however. Although we feel that the nature of the transition between the branched 
polymer phase and the inflated phase is well established, the relatively modest sys- 
tem sizes which can currently be studied make it difficult to determine v precisely, 
or even exclude the possibility that the behavior we observe in the inflated phase 
is dominated by crossover effects. Nevertheless, the fact that arguments based on 
the blob picture of de Gennes yield such a concise and consistent description of 
the observed scaling behavior makes us believe that the scenario described above 
is basically correct. In any case, a better theoretical understanding of the nature 
of the inflated phase is necessary if we are to appreciate the implications of these 
results. 
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FIGURE CAPTIONS 

1.) The average volume < V > as a function of the pressure increment p, for 
N = 47 (o), N = 127 (x), and N = 247 (0). p = Apa 3 /k B T, where the 
constant a is of the order of the interparticle spacing; for convenience we set 
a^/ksT = 1 here and in the following. 

2.) Probability distribution P(V) for N = 127, at p = 0.40 (dashed line), p = 0.50 
(full line), and p = 0.60 (dotted line). 

3.) Volume V vs. number of Monte Carlo steps per monomer, for N = 247, 
p = 0.30 and p = 0.35. In the case p = 0.30, the data are shifted by AV = —40 
in order to separate the two curves. 

4.) Scaling functions for (a) the volume, F =< V > p~ 3uj N~ Siy + , as a function 
of the scaled pressure pN 3 "/ 2 , and (b) the mean square radius of gyration, 
Sy =< R 2 g X V >~ 2 / 3 , as a function of the scaled volume < V > (N — 
No)~ 3u ' 2 , with iVo = 13. Data for both the constant-pressure ensemble with 
N = 47 (o), iV = 127 (x), and iV = 247 (D), as well as the constant volume 
ensemble with iV = 247 (o) are plotted. 

5.) Configurations of a vesicle in the two- phase region with iV = 247 monomers 
and volume V = 120, after 4, 8, 12, 16 and 20 million MCS. 

6.) The volume to area ratio jt(< V > N~ 3u ' 2 ) of an inflated vesicle. Data for 
both the constant-pressure ensemble with iV = 47 (o), iV = 127 (x), and 
iV = 247 (0), as well as the constant volume ensemble with iV = 247 (o) are 
plotted. 

7.) Scaling plot of the eigenvalue ratios 3 < Aj > / < Ai + A2 + A3 > of the 
moment of inertia tensor as a function of < V > (N — Nq)~ 3v ' 2 , with iVo = 5. 
The solid lines in (a) and (c) indicate the asymptotic scaling behavior (33) 
with the amplitudes C3 = 0.013 and c\ = —0.012, respectively. Data for both 
the constant-pressure ensemble with N — 47 (o), iV = 127 (x), and iV = 247 
(0), as well as the constant volume ensemble with iV = 247 (o) are plotted. 
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8.) Scaling plot of the asphericities Ti and 1?2 as a function of < V > JV _3l/ / 2 . 
The solid lines indicate the scaling behavior (35) with the amplitudes (cs — 
c\) = 0.025 and (03 — C2) = C3 = 0.013. Data for both the constant-pressure 
ensemble with N = 47 (o), N = 127 (x), and iV = 247 (0), as well as the 
constant volume ensemble with iV = 247 (o) are plotted. 

9.) Scaling plot of A 3 and S3 as a function of < V > N~ 3iy / 2 . Data for both the 
constant-pressure ensemble with N — 47 (o) , iV = 127 (x), and iV = 247 (0), 
as well as the constant volume ensemble with iV = 247 (o) are plotted. 
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